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1. Introduction 



Recently, there has been much theoretical interest in the possibility of using an adiabatic quantum opti- 
mization technique to solve NP-complete and NP-hard problems [H [5] . This is due to the following trick: 
suppose we have a quantum Hamiltonian Hp whose ground state corresponds to the solution of a problem 
of interest, and another Hamiltonian Hq, whose ground state is "easy" (both to find and to prepare in an 
experimental setup). Then, if we prepare a quantum system to be in the ground state of Hq, and then 
adiabatically change the Hamiltonian for a time T according to 



then if T is large enough, the quantum system will remain in the ground state for all times, and thus at 
time T will correspond to a solution of an NP problem. 

In this paper, we will focus on the first part of the problem: finding a quantum Hamiltonian, Hp, 
which can encode the ground state of a problem of interest. In particular, we will look for classical Ising 
Hamiltonians, which are extremely simple to interpret as quantum Hamiltonians by simply converting 
each classical spin variable into a qubit, and have the added advantage that they may be more adaptable 
to quantum computing hardware. In each case, there are an infinite number of possible choices of Hp 
which can be used, although many of them will be simply related to the others by, for example, slightly 
tweaking some couplings or magnetic fields. We will not ask the question of how robust the Hamiltonians 
are to slight changes in their couplings or magnetic fields, and simply content ourselves with finding a 
single Hamiltonian for each problem. 

Analogies between Ising Hamiltonians and NP problems have been frequently studied in the past 
[U |6] , although many of the NP problems previously studied are very straightforward to phrase in this 
manner. This paper can be thought of as the result of a challenge to a theoretical physicist to show how 
"all of the famous NP problems'tj [7| can be written down as Ising models. The "subtle" techniques in 
this paper are primarily of a few fiavors, which roughly correspond to the tackling the following three 
issues through a (polynomial-size) expansion of the state space: problems with inequalities as constraints 
(for example, n > 1, as opposed to n = 1), problems where the natural variables have q states (with 
q > 2), and problems which ask global questions about graphs, such as connectivity. The methods we 
use to tackle these problems with Ising models generalize very naturally to more complicated problems. 
However, it is sometimes the case that we seem to get "lucky" and find solutions which do not require 
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^No offense to anyone whose problems have been left out. 
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expansions of the state space: for example, present a solution to a tricky minimax graph problem by a 
separation of three energy scales, although this trick does not seem to generalize as well as the others. 

There has been debate about whether or not these algorithms would actually be useful for a quantum 
computer [Ull], and this paper does not settle this debate. However, it may be useful in the distant future 
to have a quantum computer solve a large NP problem without the aid of a classical computation, and so 
therefore finding quantum Hamiltonians whose ground states correspond directly to the solutions of NP 
problems may be prove useful. 

In a sense, the goal of this paper is trivial, as Ising spin glasses are known to be NP-hard [2], so it 
is natural to suspect connections with all other NP problems. This triviality is exaggerated by the fact 
that a simple NP-complete number partitioning problem is simply related to an Ising spin glass, so if 
there are polynomial time maps of all other NP-complete maps into this problem, we expect polynomial 
size Ising formulations. However, many of these direct constructions do not seem to be written down in 
a manner accessible for physicists or quantum algorithm designers, if they have been at all. Thus, the 
goal of this paper is to present these constructions in one place, in a self-contained manner, starting from 
trivial constructions and moving onwards to significantly more subtle constructions, and culminating in 
an Ising model for each of Karp's 21 NP-complete problems. 

2. Partitioning Problems 



The simplest problems to phrase as Ising models are partitioning problems: these maps are celebrated 
and well-known [8l |6]. For completeness, we review them here. Some of the techniques used in these 
models will prove useful in more nontrivial constructions. 

2.1. Number Partitioning 

Number partitioning asks the following: given a set S of numbers, which we will take to be natural 
numbers, is there a partition of this set of numbers into two disjoint subsets R and S — R, such that the 
sum of the elements in both sets is the same? This can be phrased trivially as an Ising model as follows. 
Let rii {i = 1, . . . , N = \S\) describe the numbers in set S, and let 



be an energy functional, where Sj = ±1 is an Ising spin variable. Here A > is some positive constant. 
Typically, such constants are scaled to 1 in the literature, but for simplicity we will retain them, since 
in many formulations a separation of energy scales will prove useful and retaining each scale can make it 
easier to follow. 

It is clear that if there is a solution to the Ising model with H = 0, then there is a configuration of spins 
where the sum of the rij for the +1 spins is the same for the sum of the for the —1 spins. Thus, if the 
ground state energy is H = 0, there is a solution to the number partitioning problem. Degeneracies in this 
ground state which do not correspond to s* — >• —s* correspond to distinct solutions to the partitioning 
problem. Furthermore, if the ground state has H > 0, we know that there are no solutions to the 
partitioning problem, but the ground state we do find is (one of) the best possible solutions, in the sense 
that it minimizes the mismatch. 
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2.2. Graph Partitioning 

Graph partitioning is the classic example of a map between the physics of Ising spin glasses and NP 
problems. We will phrase it in a slightly different form than the original [6], which employed a hard 
constraint on the phase space. We will want none of our formulations to do this, as this may hamper 
their ability to be used in quantum computing applications. 

Let us consider an undirected graph G = {V,E). with an even number N = \V\ of vertices. We 
ask: what is a partition of the set V into two subsets of equal size N/2 such that the number of edges 
connecting the two subsets is minimized? This is known to be an NP-complete problem. As before, we 
will place an Ising spin = ±1 on each vertex v & V on the graph, and we will let +1 and —1 denote 
the vertex being in either the + set or the — set. We solve this with an energy functional consisting of 
two components: 

H = Ha + Hb (3) 
where ^ 

is an energy which provides a penalty if the number of elements in the + set is not equal to the number 
in the — set, and 

Hb = b (5) 

{uv)€W 

is a term which provides an energy penalty B for each time that an edge connects vertices from different 
subsets. If > 0, then we wish to minimize the number of edges between the two subsets; if i? < 0, we 
will choose to maximize this number. Should we choose i? < 0, we must ensure that B is small enough 
so that it is never favorable to violate the constraint of Ha in order to minimize energy. 



2.3. Cliques 



A clique of size K in an undirected graph G = {V, E) is a subset C y of the vertices, of size \W\ = K, 
such that the subgraph {W, Ew) (where Ew is the edge set E restricted to edges between nodes in W) is 
a fully connected graph - i.e., all possible K{K — l)/2 edges in the graph are present. 

We show how the NP-complete decision problem of whether or not a clique of size K exists can be 
written as an Ising-like model. We place a spin variable St, = ±1 on each vertex v V oi the graph. In 
general, in this paper, for a spin variable Sa, we will define the binary bit variable 



So + 1 



(6) 



It will typically be more convenient to phrase the energies in terms of this variable Xa, as it will be for 



this problem. Note that any energy functional which was quadratic in will remain quadratic in 
vice versa, so we are free to use either variable. 
We will write the energy as 



and 



H = A\ K 



E 



+ A 



K{K - 1) 



(uv)eE 



(7) 



where ^ > is a positive constant. The ground state of this Hamiltonian is = if and only if a clique 
of size K exists. This is because the role of the first term is to enforce the constraint that exactly K 
vertices may be included in a trial subset W of the vertices, and the second term subsequently gives an 
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energy penalty for each possible edge between two vertices of W which is not included in E-\y- There are 
no energy penalties from the latter term precisely when the trial set W forms a clique. Given a ground 
state solution, it is of course easy to read off from the which K nodes form a clique. 

A quantum algorithm can actually be made slightly more efficient so long as the initial state can be 
carefully prepared [3]. 



3. Hard Constraint Covering and Packing Problems 



In this section, we will turn to what we call "hard constraint" problems related to covering and packing 
sets. These are essentially problems where constraints must be exactly satisfied. Many of the problems 
described below are often discussed in the literature, but again we review them here for completeness. 

3.1. Binary Integer Linear Programming 

Consider the problem, where xi, . . . , xat are binary variables, arranged into a vector x. If we ask, what is 
the largest value of c • x, for some vector c, given a constraint 



^x = b 

with A some matrix and b some vector with m components. 
We solve this as followed. Let 

H = Ha + Hb 

where 



N 



(8) 



(9) 



(10) 



and A > is a constant. The ground states of = enforce (if such a ground state exists, of course!) 
the constraint that Ax = b. Then we set 



Hb 



N 

1=1 



(11) 



with B <^ A another positive constant. The energy scale separation is required so that it is never favorable 
to violate the Ha constraint. 



3.2. Exact Cover 

The exact cover problem goes as follows: consider a set C/ = {1, . . . , n}, and subsets Vi C U such that 

U = \JVi. (12) 

i 

The question is: is there a subset of {Vi}, called R, such that the elements of R are disjoint sets, and the 
union of the elements of R is U? This problem was solved in P] but for simplicity, we repeat it here. 

We write this problem as an Ising model as follows. Let us denote a spin variable Si and binary 
variable Xi for each subset Vi. The Hamiltonian we use is H = Ha + Hb, where 

HA = Aj2 (l- E ^0 • (13) 
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Ha = precisely when every element is included exactly one time, which implies that the unions are 
disjoint. The existence of a ground state of energy H = corresponds to the existence of a solution to 
the exact cover problem. 

It is also straightforward to extend this, and find the smallest exact cover. This is done by simply 
adding a second energy scale B <^ A: 

HB = BJ2x^■ (14) 

i 

The ground state of this model will be nB, where n is the smallest number of subsets required. 
3.3. Set Packing 

Let us consider the same setup as the previous problem, but now ask a different question: what is the 
largest number of subsets Vi which are all disjoint? This is called the set packing problemU To do this, 
we use H = Ha + Hb with a separation of energy scales. We use 

Ha = A ^ XiXj, (15) 
i,j:V,nVjjtil) 

which is minimized only when all subsets are disjoint. Then, we use 

HB = -BY,^i (16) 

i 

which simply counts the number of sets we included. Choosing B < A ensures that it is never favorable 
to violate the constraint Ha (since there will always be a penalty of at least A per extra set included). 

Note that an equivalent formulation of this problem is as follows: let us consider the sets to be encoded 
in an undirected graph G = (y,E), where each set Vi corresponds to a vertex i € V. An edge ij € E 
exists when Vi n Vj is nonempty. It is straightforward to see that if we replace 

Ha = AY, XiXj (17) 
ijeE 

that the question of what is the maximal number of vertices which may be "colored" {xi = 1) such that no 
two colored vertices are connected by an edge, is exactly equivalent to the set packing problem described 
above. 



3.4. Vertex Cover 

Given an undirected graph G = {V,E), what is the smallest number of vertices that can be "colored" 
such that every edge is incident to a colored vertex? Let x^ be a binary variable on each vertex, which is 
1 if it is colored, and if it is not colored. Our Hamiltonian will he H = Ha + Hb- The constraint that 
every edge has at least colored vertex is encoded in Ha'- 

Ha = A^{1-Xu){1-x,). (18) 

uv&E 

Then, we want to minimize the number of colored vertices with Hb- 

Hb = bY,Xv (19) 

V 

Choose NB < A, where N = \V\. Then the ground state of this Ising model corresponds to the config- 
uration of colored vertices such that each edge connects to a colored vertex, and the fewest number of 
vertices are colored. 

^Often times, this is also referred to as the maximaUy independent set (MIS) problem. 
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3.5. Satisfiability 

Satisfiability is one of the most famous NP-complete problems. Every satisfiability problem can be written 
as a so-called 3SAT problem in conjunctive normal form (and this algorithm takes only polynomial 
steps/time) and so we will focus for simplicity on this case. In this case, we ask whether 

^ = Ci A C72 • • • A C7„ (20) 

can take on the value of true - i.e., every Cj for 1 < -i < is true, where the form of each Ci is: 

Ci = y^i V V (21) 

Here y^j, and are selected from another set of Boolean variables: xi, . . . ,xj^ ,xi, . . . ,xn- 

There is a well-known reduction of 3SAT to the set packing problem [3] which we reproduce here, for 
completeness. Consider solving the set packing problem on the following graph G with 3m nodes, which 
we construct as follows. For each clause Ci, we add 3 nodes to the graph, and connect each node to the 
other 3. After this step, if there is a yi and y2 such that yi = Xj and y2 = xj, then we also add an edge 
between these two nodes. Solving the set packing problem on this G, and asking whether the solution has 
exactly m nodes, is equivalent to solving the 3SAT problem. This can be seen as follows: if a solution to 
the 3SAT problem exists, only one element of each clause needs to be true - if more are true, that is also 
acceptable, but we must have that one is true, so let us choose to color the vertex corresponding to the 
variable which is true. However, we may also not choose to have both xi be true and xi be true, so we 
are required to connect all such points with an edge. 



3.6. Minimax Graph Problem 

The minimax graph problem proceeds as follows: let G = {V, E) denote an undirected graph, and let 
C C E' be a proposed "coloring". The constraints on C are as follows: for each edge in C, let us color the 
two vertices it connects. We will then demand that: no two edges in C share a vertex, and that there are 
no two uncolored vertices connected by an edge. This is minimal in that we cannot add any more edges 
to C (coloring any appropriate vertices) without violating the first constraint, and maximal in the sense 
that the trivial empty set solution is not allowed - we must include all edges between uncolored vertices. 

Note that, from this point on in this paper, we have not found the Ising formulations of this paper in 
the literature. 

We will use the spins on the graph to model whether or not an edge is colored. Let us use the binary 
variable Xe to denote whether or not an edge is colored. Our method consists of the following trick. Our 
energy will consist of a series of three energy functionals, each of which is chosen to be "much" smaller 
than the previous one: 

H = Ha + Hb + He. (22) 

The first and largest term, H^, will impose the constraint that no vertex has two colored edges. This can 
be done by setting 

Ha = AY^ Xe.Xe,. (23) 

^ {ei,e2}Cdv 

Here ^4 > is a positive energy, and dv corresponds to the subset of E of edges which connect to v. Thus 
the ground states consist of Ha = 0; if Ha > 0, it is because there is a vertex where two of its edges are 
colored. 

If we choose A to be large enough, then we also can define, for states with Ha = 0, the variable 



Vv 



1 V has a colored edge 
V has no colored edges 



E (24) 
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We stress that this definition is only vahd for states witli = 0, since in these states each vertex has 
either or 1 colored edges. We then define the energy Hb, such that solutions to the minimax coloring 
problem also have Hb = 0. Since we have already constrained the number of colored edges per vertex, 
we choose Hb to raise the energy of all solutions where there exists a possible edge which can be colored, 
yet still not violate the coloring condition, out of the ground state. To do this, we can sum over all edges 
in the graph, and check whether or not the edge connects two vertices, neither of which are colored: 

Hb = b J2 {i-yu){i-yv). (25) 

e={uv) 

Note that since, can be negative, wc must choose i? > to be small enough. This will ensure that a 
ground state will have energy Ha + Hb = 0, and correspond precisely to Ha = Hb = 0: i.e., states which 
do not violate the minimax constraints. 

Now, given the states where Ha = Hb = 0, we now want the ground state to be the state where the 
fewest number of edges are colored. To do this, we simply let 

Hc = Cj2xe (26) 

e 

count the number of colored edges. Here C is an energy scale chosen to be small enough such that it 
is never energetically favorable to violate the constraints imposed by either the Ha or Hb terms: for 
example, CN <2B. The term with the smallest He has the smallest number of edges, and is clearly the 
solution to the minimax problem. Each ground state of this spin model is equivalent to a solution of the 
minimax problem. 



4. Soft Constraint Packing and Covering Problems 



We now turn to NP problems whose formulations as Ising models are substantially more subtle. This is, 
in every case, due to the fact that the constraints on the system are "soft" constraints which are only 
inequalities. For most problems of interest, these soft constraints can be re-written as "hard constraints" 
by an expansion of the number of spins. 

4.1. Set Cover 

Consider a set U = {!,..., n}, with sets Va QU (a = 1, . . . , A'') such that 

N 

U=[^Va. (27) 

a=l 

The set covering problem is to find the smallest possible number of V^s, such that the union of them is 
equal to U . 

Our solution consists of the following. Let us denote Xa to be a binary variable which is 1 if set a is 
included, and if set a is not included. Let us then denote Xi^m to be a binary variable which is 1 if the 
number of V^s which include element i is m > 1, and otherwise. Set H = Ha + Hb- Our first energy 
imposes the constraints that exactly one Xi^m must be 1, since each clement of U must be included a fixed 
number of times, and that the number of times that we claimed i was included is in fact equal to the 
number of we have included, with i as an element: 

n / N \2 n / N \2 

HA = AY^\l-Y^XiA +A^{Y^mxi^m- ^ xA . (28) 

1=1 V m=l / i=l \m=l a:i€V"a / 
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Finally, we minimize over the number of V^s included: 

HB=BY,Xa, (29) 
a 

with B <^ A chosen so that the A constraint is never violated. 
4.2. Knapsack with Integer Weights 

The knapsack problem is the following NP-complete problem: given a list of N objects, labeled by indices 
a, with the weight of each object given by Wa, and a knapsack which can only carry weight W, what is 
the heaviest collection of objects which we can add to the knapsack such that their combined weight is 
smaller than Wl 

Let Xa be a binary variable denoting if we included a, and x„ for 1 < n < W denote a binary variable 
which is 1 if the final weight of the knapsack is n, and otherwise. Our solution consists of letting 
H = Ha + Hb, with 

1 - ^ Xn J + ^ X] nXn - ^ WaXa 1 (30) 
n=l / \n=l a / 

which enforces that the weight can only take on one value and that the weight of the objects in the 
knapsack equals the value we claimed it did, and finally 

w 

HB = -Bj2nxn, (31) 

n=l 

with B > and B <^ A, chosen such that it is never favorable to violate the constraints. 

There is a trick which can be used to dramatically reduce the number of extra Xi spins which must 
be added. For simplicity, let us focus on the case W = 2^, for M some positive integer, where the trick 
works most efficiently. In this case, we only need xi, . . . ,xm instead of Xi,. . . ,xw- It is easy to check 
that 

(M 
^2^-'xn-^WaxA -B^2^-^Xn (32) 
n=l a / n=l 

solves the exact same knapsack problem. This trick can be used on many of the problems below as well, 
but it makes the physical intuition a bit less clear, so we will not use it in our explicit formulations. 



5. Coloring Problems 



We now turn to coloring problems. Naively, coloring problems are often best phrased as Potts models, but 
these classical Potts models can be converted to classical Ising models with an expansion of the number 
of spins. This simple trick forms the basis for our solutions to this class of problems. 

5.1. Graph Coloring 

Given an undirected graph G = {V, E), and a set of n colors, is it possible to color each vertex in the graph 
with a specific color, such that no edge connects two vertices of the same color? This is a generalization 
of the classic problem of how many colors are needed to color a map, such that every two countries which 
share a border have a different color, and is called the graph coloring problem. 
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Our solution consists of the following: we denote Xy^i to be a binary variable which is 1 if vertex v is 
colored with color i, and otherwise. The energy is 



(33) 



i=l 



The first term enforces the constraint that each vertex has exactly one color, and provides an energy 
penalty each time this is violated, and the second term gives an energy penalty each time an edge 
connects two vertices of the same color. If there is a ground state of this model with H = 0, then there 
is a solution to the coloring problem on this graph with n colors. We can also read off the color of each 
node (in one such coloring scheme) by looking at which xs are 1. Note that the number of spins can be 
slightly reduced, since there is a permutation symmetry among colorings, by choosing a specific node in 
the graph to have the color 1, and one of its neighbors to have the color 2, for example. 

5.2. Clique Cover 

The clique cover problem, for an undirected graph G = {V,E), is the following: given n colors, is there a 
coloring of the graph sTich that if we assign every vertex in the graph exactly one color, then the subgraph 
for any color is a clique? 

Our solution is very similar to the graph coloring problem. Again, we employ the same binary variables 
as for graph coloring, and use Hamiltonian 



H 



V V / V {uv)eW 



(34) 



The first term enforces the constraint that each vertex has exactly one color by giving an energy penalty 
each time this constraint is violated. In the second term, since the sum over v of x^^i counts the number 
of nodes with color i, the first sum counts highest possible number of edges that could exist with color i. 
The second term then checks if, in fact, this number of edges does in fact exist. Thus = if and only 
if the clique cover problem is solved by the given coloring. If a ground state exists with H = 0, there is a 
solution to the clique covering problem. 



5.3. Job Sequencing with Integer Lengths 

The job sequencing problem is as follows: given a list of N jobs for, say, a computer cluster, and job i has 
length Li, how can each job be assigned to a computer in the cluster such that, if the length of cluster a, 
and its set of jobs is Va, then 

We assume that G N. 

To do this, we will begin by demanding that without loss of generality, Mi > for any a. Introduce 

the variables Xj^„ which are 1 if job i is added to computer a, and otherwise, and the variables yn,a for 
a 7^ 1 and n > 0, which is 1 if the difference M\ — = n. Then the Hamiltonian 

N / \2m/ \2 

Ha = A^il -^Xi^aj + A^i^nyn,a + ^Li{xi^a- Xi,l)\ (36) 
i=l \ a / Q=l \ n i / 
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encodes that each job can be given to exactly one computer, and that no computer can have a longer 
total length than computer 1. To find the minimal maximal length Mi, we just use 

HB = Bj2LiXi,i. (37) 

i 

6. Hamiltonian Cycles 



In this section, we describe the solution to the (undirected or directed) Hamiltonian cycles problem, 
and subsequently the traveling salesman problem, which for the Ising spin glass formulation, is a trivial 
extension. 

6.1. Hamiltonian Cycles and Paths 

Let G = {V,E), and = \V\. The graph can either be directed or undirected, and the solution will not 
change. The Hamiltonian path problem is as follows: starting at some node in the graph, can one travel 
along an edge, visiting other nodes in the graph, such that one can reach every single node in the graph 
without ever returning to the same node twice? The Hamiltonian cycles problem asks that, in addition, 
the traveler can return to the starting point from the last node he visits. 

Without loss of generality, let us label the vertices 1, . . . , AT, and take the edge set (uv) to be directed 
— i.e., the order uv matters. It is trivial to extend to undirected graphs, by just considering a directed 
graph with (vu) added to the edge set whenever (uv) is added to the edge set. Our solution will use N'^ 
spins: binary bit variables Xy^i, where v represents the vertex and i represents its order in a prospective 
cycle. Our energy will have three components. The first two things we require is to enforce that every 
vertex can only appear once in a cycle, and that there must be a j^^ node in the cycle for each j. Finally, 
for the nodes in our prospective ordering, if x^j and x^j+i are both 1, then (uv) € E. Note that N + 1 
should be read as 1, in the expressions below, if we are solving the cycles problem. These are encoded in 
the Hamiltonian: 

n / N \^ n / N N 

f=l \ J=l / j=l V V=l ) (u„)^£j=l 

A > is a constant. It is clear that a ground state of this system has H = only if we have an ordering 
of vertices where each vertex is only included once, and adjacent vertices in the cycle have edges on the 
graph - i.e., we have a Hamiltonian cycle. 

To solve Hamiltonian path, instead of Hamiltonian cycle, all we have to do is restrict the sum over 
j above from 1 to iV — 1. Thus it does not care about whether or not the first and last nodes are also 
connected. 

We note that it is straightforward to imagine slightly reducing the size of the state space for the 
Hamiltonian cycles problem as follows: it is clear that node 1 must always be included in a Hamiltonian 
cycle, and without loss of generality we can set xi^i = 1: this just means that the overall ordering of the 
cycle is chosen so that node 1 comes first. This reduces the number of spins to (A'" — 1)^. 

6.2. Traveling Salesman 

The traveling salesman problem for a graph G = {V,E), where each edge uv in the graph has a weight 
Wuv associated to it, is to find the Hamiltonian cycle such that the sum of the weights of each edge in 
the cycle is minimized. Typically, the traveling salesman problem assumes a complete graph, but we have 
the technology developed to solve it on a more arbitrary graph. 
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To solve this problem, we use H = Ha + Hb, with Ha the Hamiltonian given for the directed (or 
undirected, if the graph is undirected for traveling salesman) Hamiltonian cycles problem. We then simply 
add 

N 

Hb = B Wu,Yl 

(39) 

{uv)&E j = l 

with B small enough that it is never favorable to violate the constraints of Ha- If the traveling salesman 
does not have to return to his starting position, we can restrict the sum over j from 1 to A?^ — 1, as before. 



7. Tree Problems 



The most subtle NP problems to solve with Ising models are problems which require finding connected 
tree subgraphs of larger graphs. The key point which makes the solution difficult is the requirement 
of connectivity, and the solution relies on similar tricks, therefore, to the Hamiltonian cycles problem 
discussed earlier. 

7.1. Minimal Spanning Tree with a Maximal Degree Constraint 

The minimal spanning tree is the following: given an undirected graph G = {V^E), where each edge 
{uv) € E comes with a cost Cuv, what is the tree T C G, which contains all vertices, such that the cost of 
T, defined as 

c{T)= ^ c„^, (40) 

is minimized? Without loss of generality, we take c^v > in this subsection (a large positive constant can 
always be added to each c^v ensure that the smallest value of c^^ is strictly positive, since the number of 
edges of a tree is always — 1, if the number of vertices in N). We will also add a degree constraint, 
that each degree in T be < A, which makes the problem NP-complete. 

To solve this problem, we place a binary variable ye on each edge to determine whether or not that 
edge is included in T: 

~ \ otherwise ' ^ ' 

We also place a large number of binary variables x^^i on each vertex, and Xuv,i 7^ Xvu,i on edge (uv): the 
number i = 0, 1, ... , N/2 will be used to keep track of the depth a node in the tree, and if Xuv = 1, it 
means that u is closer to the root than v, and if x^,u = 1 it means that v is closer to the root, and we 
use another variable z^^i (i = 1, ... A) to count the number of degrees of each node. We now use energy 
H = Ha + Hb, where the terms in Ha are used to impose the constraints that: there is exactly one root 
to the tree, each vertex has a depth, each bond has a depth, and its two vertices must be at different 
heights, the tree is connected (i.e., exactly one edge to a non-root vertex comes from a vertex at lower 
depth), each node can have at most A edges, and each edge at depth i points between a node at depth 
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i — 1 and i, respectively: 

Ha = Ail -'^Xyfi] + I 1 -^a;„,i I + ^ ^ I y„„ - ^(a;„^,i + a;^„,i) 1 

\ V / V \ i / uveE \ i / 

N/2 / \ ^ N/2 / A 

+ ^ X! X! I ^'"''^ ~ XI I ^ XI XI I XI ^^'"'^ ~ XI XI 

V i=l y {uv)eE j V 1=1 yj=l uv,vueE i 

N/2 

UV,VUS:E 4 = 1 

The ground states with //^ = arc trees which include every vertex. Note that when we sum over uv G E, 
we are not counting vu as distinct, but when we notate uv, vu € E, we do treat uv ^ vu. We then add 

N/2 

(43) 

uv,vu&E i=\ 

The minimum of E will find the minimal spanning tree, subject to the degree constraint. 

7.2. Steiner Trees 

The Steiner tree problem is somewhat similar to the problem above: given our costs Cuv, we want to find 
a minimal spanning tree for a subset U <ZV o{ the vertices, with no degree constraints. To do this, we 
use the same Hamiltonian as for the minimal spanning tree, except we add binary variables for v ^ U 
which determine whether or not a node U is included, and use the Hamiltonian H = Ha + Hb, where 
Ha enforces constraints similarly to in the previous case: 

Ha = All - ^ X^,0 I + ^ XI ( ~ XI ^'''^ ) [Vuv- E{Xuv,i + Xyu,i) J 

\ V / V \ i / uv&E \ i / 

N/2 / \ ^ N/2 

~^ ^ ^ ^ ^ ^ I Xy^i ^ ^ Xuv,i J ~l~ A ^ ^ ^ ^ Xuv,i{'2: Xu,i—1 Xy^i) (44) 

V i=l y {uv)eE j uVjVueE i=l 

We then use Hb from the previous model to determine the minimum weight tree. 

7.3. Directed Feedback Vertex Set 

A feedback vertex set for a directed graph G = (V, E) is a subset F d X such that the subgraph 
{X — F,d{X — F)) is acyclic. We will refer to F as the feedback set. Solving a decision problem for 
whether or not a feedback set exists for \F\ < k is NP-complete. We solve the optimization problem of 
finding the smallest size of the feedback set first for a directed graph - the extension to an undirected 
graph will be a bit more involved. 

Before solving this problem, it will help to prove two lemmas. The first lemma is quite simple: there 
exists a node in a directed acyclic graph which is not the end point of any edges. For suppose that for 
each vertex, there was an edge that ends on that vertex. Then pick an arbitrary vertex, pick any edge 
ending on that vertex, and follow that edge in reverse to the starting vertex. Repeat this process more 
than N times, and a simple counting argTimcnt implies that wc must have visited the same node more 
than once, at least once. Thus, we have traversed a cycle in reverse, which contradicts our assumption. 
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The second lemma is as follows: a directed graph G = (V, E) is acyclic if and only if there is a height 
function h : V —i- N such that if uv G E, h{u) < h{v): i.e., every edge points from a node at lower 
height to one at higher height. That height function existence implies acyclic is easiest to prove using the 
contrapositive: suppose that a graph is cyclic. Then on a cycle of edges, we have 



is a contradiction. To prove that an acyclic graph has a height function, we construct one recursively. 
Using our first lemma, we know that there exists a vertex u with only outgoing edges, so let us call 
h{u) = 1. For any other vertex, we will call the height of that vertex h{v) = 1 + h'{v), where h'{v) is 
found by repeating this process on the graph with node u removed (which must also be acyclic). It is 
clear this process will terminate and assign exactly one node height i for each integer 1 < i < 

We can now exploit this lemma to write down an Ising spin formulation of this problem. We place a 
binary variable on each vertex, which is if u is part of the feedback set, and 1 otherwise. We then 
place a binary variable Xy^i on each vertex, which is 1 if vertex v is at height i. So far the heights i are 
arbitrary, and the requirement that a height function be valid will be imposed by the energy. The energy 
functional we use is if = Ha + Hb where 



The first term ensures that if a vertex is not part of the feedback set, it has a well-defined height; the 
second term ensures that an edge only connects a node with lower height to a node at higher height. We 
then find the smallest possible feedback set by adding 



7.4. Undirected Feedback Vertex Set 

The extension to undirected graphs requires a bit more care. In this case, we have to be careful because 
there is no a priori distinction on whether the height of one end of an edge is smaller or larger than the 
other - this makes the problem much more involved, at first sight. Furthermore, it is not true that a 
directed acyclic graph is acyclic if the orientation of edges is ignored. However, for an undirected graph, 
we also know that a feedback vertex set must reduce the graph to trees, although there is no longer 
a requirement that these trees are connected (this is called a forest). With this in mind, we find that 
the problem is actually extremely similar to minimal spanning tree, but without degree constraints or 
connectivity constraints. The new subtlety, however, is that we cannot remove edges. 

To solve this problem, we do the following: introduce a binary variable x,i,^i, which is 1 if v is a vertex 
in any tree at depth i, and otherwise. However, to account for the fact that we may remove vertices, 
we will allow for x^^_i = 1 if v is part of the feedback vertex set, and otherwise. We do a similar thing 
for edges: we consider Xuv,i 7^ Xvu,i to be defined as before when i > 0. We also define the variables 
Xuv,-i 7^ Xyu,-i^ which we take to be 1 when the ending node of the "directed" edge is in the feedback 
vertex set. Now, we can write down a very similar energy to the minimal spanning tree: 



< ^[/l('Ui+l) - h{ui)] = h{ui) - h{un) + h{un) - h{un-i) H h{ui) = 



(45) 




(46) 




(47) 



V 




(48) 
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The changes are as fonows: we no longer constrain only 1 node to be the root, or constrain the degree of a 
vertex - however, we have to add a new term to ensure that edges are only ignored in the tree constraint 
if they point to a node in the feedback set. We then add 



Hb = B^x,,.i (49) 



with B <^ A chosen so that the A constraints are never violated. This counts the number of nodes in the 
feedback set, so thus H is minimized when Hb is smallest - i.e., we have to remove the fewest number of 
nodes. 

7.5. Feedback Edge Set 

For a directed graph, the feedback edge set problem is to find the smallest set of edges F C E such that 
(V, E — F) is a directed acyclic graph. It is known to be NP-hardU Our solution will be somewhat similar 
to the directed feedback vertex set. We place a binary variable yuv on each edge, which is 1 if uv ^ F, and 
define Xuv,i to be 1 if uv if both Uuv = 1 and the height of node u is i. We also add a binary variable x^^i, 
as for the feedback vertex set. Our constraint energy must then enforce that: each vertex and included 
edge has a well-defined height, and that each edge points from a lower height to a higher height: 

ha = a'^Ii Xy^i + ^ X] h^™ ~ X] + ^ ~ ~ X^ ^^'-j' ■ 

f \ i / uv<=E \ i / uv i \ j>i J 

We then use 

HB = BY,i'^- yuv) (51) 

uv£E 

to count the number of edges in i*" - it is minimized when this number is smallest. 



8. Graph Isomorphisms 



The question of whether two graphs Gi = {Vi,Ei) and G2 = {V2,E2) are isomorphic is believed to be 
hard, but its classification into a complexity class is still a mystery. Since it is a hard problem, let us 
nonetheless describe an Ising formulation for it. An isomorphism is only possible if \Vi\ = IV2I ^ N, so 
we will restrict ourselves to this case, and without loss of generality, we label the vertices of Gi with 
l,...,iV. 

We write this as an Ising model as follows. Let us describe a proposed isomorphism through binary 
variables x^j.i which is 1 if vertex v in G2 gets mapped to vertex i in Gi. The energy 

ensures that this map is bijective. We then use an energy 

Hb = B ^ 1 - ^ Xu,iXyj \ + B ^ ^ Xu,iXvj (53) 

ijeEi y UV&E2 j ij^Ex UV&E2 

to penalize a bad mapping: i.e. an edge that is not in Gi is in G2, or an edge that is in Gi is not in G2. 
If the ground state of this Hamiltonian corresponds to = 0, there is an isomorphism. 



^It is in P if the graph is undirected however. 
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In this paper, we have presented classical Ising formulations for a wide variety of famous NP problems. 
Although most of the new constructions are quite cumbersome and unlikely to be useful for a long time, 
if ever, the tricks introduced may be able to be refined and made efficient enough to become useful. Much 
of this paper can be thought of as an amusing exercise to simply enumerate many subtle constructions, 
although it may also turn out that some of the new constructions are quite useful for quantum algorithm 
designers as well. We stress that techniques involving separations of multiple scales are quite likely the 
most efficient for solving nontrivial problems, and this may be a useful direction to explore further. 

We note that many of the Ising formulations of tree problems in this paper required substantial 
expansions of the state space and perhaps there are ways to avoid this with more clever Hamiltonians. I 
am unsure if it is an open question to determine the smallest possible scaling of state space size which 
can encode arbitrary instances of many NP problems, and this may be a worthwhile future direction. 

The fact that there is an Ising model representation for a problem by no means ensure it is NP- 
complete or NP-hard. As a reminder that sometimes the Ising approach may be quite inconvenient, let 
us consider the following example. Consider the simple problem of finding the largest integer in a list 
ni, . . . , njv. Introducing binary variables Xj for i = 1, . . . , N , the Ising model 



for A > Bmax{ni) solves this problem. In fact, this problem looks somewhat like an instance of the 
random field Ising model on a complete graph, and yet this has a very simple 0{N) classical algorithm. 
Thus, the results of this paper should be taken with a grain of salt - it is quite possible there are far more 
efficient algorithms, even for a quantum computer. 
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